Plasmon tunability in metallodielectric metamaterials 
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The dielectric properties of metamaterials consisting of periodically arranged metallic nanoparti- 
cles of spherical shape are calculated by rigorously solving Maxwell's equations. Effective dielectric 
functions are obtained by comparing the reflectivity of planar surfaces limiting these materials with 
Fresnel's formulas for equivalent homogeneous media, showing mixing and splitting of individual- 
particle modes due to inter-particle interaction. Detailed results for simple cubic and fee crystals 
of aluminum spheres in vacuum, silver spheres in vacuum, and silver spheres in a silicon matrix are 
presented. The filling fraction of the metal / is shown to determine the position of the plasmon 
modes of these metamaterials. Significant deviations are observed with respect to Maxwell-Garnett 
effective medium theory for large /, and multiple plasmons are predicted to exist in contrast to 
Maxwell-Garnett theory. 
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I. INTRODUCTION 

Artificial materials with tailored optical response func- 
tions (metamaterials) constitute a rich field of research 
due to their applicability to facilitate the design of highly- 
demanding optical devices for computing and communi- 
cations technologies. That is the case of left-handed ma- 
terials, which posses negative index of refraction^ and 
permit designing lenses of resolution below the diffrac- 
tion limitjS Furthermore, newly available nanostructured 
materials^i^:^ pose interesting possibilities that deserve 
detailed exploration to search for extraordinary optical 
properties. 

Composite materials with components of dimensions 
much smaller than the light wavelength can be assimi- 
lated in general to equivalent homogeneous media with 
macroscopic optical properties such as a dielectric func- 
tion and a magnetic permeability. The problem of finding 
the effective optical response of arrays of homogeneous 
spheres was already discussed by Maxwell.^ In fact, ef- 
fective medium theories that yield simple analytical ex- 
pressions have been available for nearly a century, like 
in the case of Maxwell-Garnett 's theory (MGT))^ which 
is particularly suited to describe spherical inclusions and 
rather accurate at small concentrations of the latter /, 
or Bruggeman's theory;^ less accurate but intended to 
work for arbitrary values of /. Extensions and correc- 
tions to these theories have been reported over the last 
decades for periodioia.ii^'i^^'i^^'ii^ia^^^ and 
disordere d^^d^i^'^1^'^'^^1^^'^'' composites, all of them shar- 
ing in common the assumption of local response. 

For metallic components plasmons are observed as col- 
lective oscillations of the valence electrons. In bulk met- 
als, they are signalled by the vanishing of bulk dielectric 
function, e = 0. For an ideal metallic behavior as that 
described by Drude's formula 
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with 77 ^ 0+, this occurs at the bulk plasmon frequency 
uj ^ LUp. The plasmon frequency is affected by the 
shape of the boundary of the metal, and for instance 
surface plasmons occur with frequency ujg — i^p/V^ at 
a planar surface, whereas modes of frequencies lui = 
ujpy^l/{2l + I) are observed in spherical particles in the 
long- wavelength limit, corresponding to multipole oscilla- 
tions of orders Z = 1, 2, . . . , as described in the early work 
of Mie^ Multipolar plasmons can be also found in non- 
spherical metallic particles like nanorods,^^ nanorings,^^ 
and nanoshells.'^^ Collective modes in these particles and 
in particle pairs have been explained intuitively by in- 
voking plasmon chemistry concepts following the gen- 
eral procedure outlined in Ref. l32l For layers of par- 
ticles on a substrate, new plasmons are observed as a 
result of their mutual interactioni^ In 3D periodic ar- 
rangements of spheres like the ones discussed below, In- 
glesfield et al?^ have reported a comprehensive study of 
plasmon bands and have found that wide bands do ex- 
ist and that significant deviations from MGT occur in 
near-touching metal spheres. This connects to previous 
results on plasmon-like behavior in wired structures at 
THz frequencies 1^ 

In this work we concentrate on metamaterials formed 
by periodic arrays of metallic spherical nanoparticles con- 
fined within a host dielectric medium. We calculate 
their effective dielectric function from a rigorous solution 
of Maxwell's equations, therefore including all multipole 
corrections in an exact fashion, as a function of metallic 
filling fraction /. Tunable plasmons like the ones under 
study here can find application to light sensors, absorbers 
and emitters of light, and functional lasers The rele- 
vance of this study is illustrated by the several methods 
that have been developed for synthesizing nanocompos- 
ites, making use of nanoparticles self-assembly2i^i22:^ to 
prepare ordered layers of metallic particles in dielectric 
hosts, including fee arrangements of Au particles in Si2^ 
and other complex particle shapes 1^1^ 

This paper shows how the plasma frequencies can be 
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engineered in metallodielectric composites consisting of 
spherical metallic nanoparticles embedded in a dielectric 
host. We consider in particular simple cubic (sc) and 
fee arrays of aluminium nanoparticles in vacuum, as well 
as fee arrays of silver nanoparticles surrounded by either 
vacuum or silicon. Modes other than dipolar are not 
resolved by external light in isolated nanoparticles, which 
leads to the prediction of a single plasmon branch as a 
function of metallic filling fraction when MGT is used. 
However, our rigorous solution of Maxwell's equations 
reveals several plasmon branches as a result of mixing 
and splitting of individual-particle modes of multipolar 
nature. 

II. PLASMON CHEMISTRY 

The basic ingredients of the composites under consid- 
eration are metallic nanoparticles. A isolated metallic 
nanoparticle can exhibit dipolar excitations with triple 
degeneracy (along the three orthogonal directions of 
space). When the dielectric function of the metal is 
described by Eq. dT]), the frequency of these modes is 
uji = u!p/^/3^ However, the interaction between neigh- 
boring nanoparticles can produce mixing and splitting 
of these dipolar modes, in a way similar to the mixing 
and splitting of p orbitals in molecular binding, which 
suggests the term of plasmon chemistry. The mode fre- 
quencies and their corresponding polarization patterns 
for two spheres of radius R interacting at a certain dis- 
tance d are represented in Fig. [Jl where the two in- 
termediate frequency states are doubly degenerate. Of 
these modes, on ly the lowest lying one and the third one 
{ojiy^l + / d?) display a net dipole moment, and there- 
fore, external light with these two frequencies should cou- 
ple strongly to this structure, whereas the other modes 
are either symmetry forbidden, depending on the light in- 
cidence conditions, or contribute very little to light scat- 
tering (silent modes). 

These conclusions are corroborated by the scatter- 
ing cross section represented in Fig. [2] for two small Al 
spheres as a function of light frequency and distance be- 
tween the sphere centers d. Two different polarization 
directions have been considered, as shown in the insets. 
Results similar to those of Figs. [1] and have been re- 
ported in Ref. |3^. 

The cross section is obtained by using the multiple elas- 
tic scattering of multipole expansions method (MESME), 
designed to solve Maxwell's equations in the presence 
of large numbers of particles in arbitrary positions4^ In 
it, the electric field is expressed in frequency space to in 
terms of magnetic and electric scalar functions 'ipf^ cind 
V'f', respectively, as^ 

E = ^[L,V^f -AvxL,^f], 

i 

where fc = w/c, = — i(r — r^) x V is the orbital angular- 
momentum operator, the sum is extended over particle 
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FIG. 1: Plasmon chemistry in interacting nanoparticles. The 
non-retarded plasmon energies of two identical metallic par- 
ticles are shown along with schematic representations of the 
corresponding oscillation modes. The particles have radius R, 
their center-to-center distance is d, and the metal is described 
by a Drude dielectric function of bulk-plasma frequency ojp 
[Eq. (UJI 
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FIG. 2: Total scattering cross section cr of a system formed 
by two metallic spheres as a function of incident photon en- 
ergy Lu and center-to-center distance d between two aluminum 
spheres described by the Drude dielectric function of Eq. ^ 
with LUp = 15 eV and rj = 0.6 eV. Two different polarization 
directions of the external electric field have been considered 
in (a) and (b), as shown in the insets, a is normalized to the 
projected area of the two spheres, 2txR? . The solid curves 
represent the modes given by the expressions of Fig. [T] 

positions r^, and M and E refer to magnetic and elec- 
tric polarization, respectively. The electromagnetic field 
created by some external source is expanded in terms of 
spherical waves as 

(r) = E i'^' - 1)^'" (?^) ' 

hn 

where ji is a spherical Bessel function, Yim is a spherical 
harmonic, and a polarization index (either magnetic or 
electric) is understood. In the absence of multiple scat- 
tering, the induced part of the scalar functions ■0™'^ = ■01'' 
results from single scattering (ss) of -0™* by each object: 

(r) = Y^i'h^+\k\v-v,\)YU^^)^ti^. (2) 

Ini 

for r outside the particle centered at r^. Here, h'f^^ is 
a spherical Hankel function. The relation between 0^™' 
and Tpl^ is provided by the scattering matrix ij, implic- 
itly defined by ■0f'' = ti'tjjf^*' , where ■0^''^'"^*-' is a vector 



of components V'r/jn'*^ ■ ^he matrix U describes the full 
scattering properties of object r^. In a cluster of several 
particles, V™'^ takes the same form as Eq. except 
that it is made up of (single scattering of the exter- 
nal field at r^) plus the result of the free propagation of 
^ind fj.Q]-^ each object rj ^ r^, followed by scattering at 
Ti (self-consistent multiple scattering). That is;^ 
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where the operator Hij describes the noted propagation. 
In this work we consider homogeneous spheres, so that ti 
becomes diagonal and given by analytical expressions^^ 
Furthermore, we have found convergent results, using a 
finite number of multipoles with I < 8. 

For a bi-sphere system, Fig. [DJb) demonstrates that 
only the lowest-frequency mode is excited for incident- 
light polarization along the inter-particle direction, 
whereas Fig. [D^a) shows that only the third hybridized 
dipolar mode contributes for polarization perpendicular 
to the inter-particle direction, as expected from the above 
discussion. The plasmon modes derived from the simple 
dipole-mode hybridization model of Fig. [1] follow the re- 
gions of large cross section only for relatively large val- 
ues of d (solid curves in Fig. [2]). Moreover, new reso- 
nances show up for parallel polarization [Fig. [D(b)] when 
the spheres are close together, as a result of stronger in- 
teraction among multipoles of the spheres that produce 
hybridized modes with net dipole moment along the ex- 
ternal field direction. 

In a structure formed by a distribution of small 
nanoparticles, one can use the Clausius-Mossotti 
formula''^ that links the polarizability of the particles a 
to the effective dielectric function of such material Coff: 



£cff — eh 
Ceff + 2eh 



An a 
3eh V ' 
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where v is the average volume per particle and eh is the 
dielectric function of the host material. This formula 
does not take into consideration non-dipolar interactions, 
and the effect of local order is neglected. However, it 
works quite well for small particle concentrations. The 
polarizability of metallic nanoparticles of radius R can 
be obtained as a function of the dielectric function of the 
metal en 
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The combination of Eqs. (|3|) and ^ is equivalent to 
MGT4 This theory predicts a single plasmon for the 
composite medium that obeys the relation 
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as a function of metal filling fraction / = AttR^/Sv. Eq. 
^ has been represented in Figs. [3] and d] by solid curves. 
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FIG. 3: Contour plot of Im{ — |e(,ft|/ecff} as a function of in- 
cident photon energy uj and filling fraction of the metal for 
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FIG. 4: Contour plot of Im{ — |eefr|/ecfi } as a function of inci- 
dent photon energy ui and filling fraction of the metal for silver 
spheres in an fee configuration, surrounded by (a) vacuum 
or (b) silicon. Brighter regions correspond to higher values 
of Im{ — |eeff|/ecff }• The solid curves correspond to Maxwell- 
Garnett theory as given by Eq. 



When the spheres are relatively close to each other, 
multipolar terms become relevant in the effective re- 
sponse, which are not accounted for by Eq. Here, we 
have included them by rigorously solving Maxwell's equa- 
tions using the MESME formalism sketched above and 
summarized in Eq. ([3]), as applied to periodic composites 
consisting of a large number of layers of nanoparticles. 
First, reflectance and transmittance matrices of incom- 
ing and scattered plane waves are obtained for each layer 
using a MESME-like procedure. Then, the interaction 
among layers is described by combining their reflectance 
and transmission matrices. More details about this so- 
called layer-KKR method, along with an efficient imple- 
mentation, can be found in the excellent work of Stefanou 
et o/4^ We have followed their approach and have used 
less than 151 plane waves to obtain convergence in the 
examples that follow, with the largest number of waves 
needed when the spheres are nearly close-packet, as ex- 
pected from previously reported divergence effects.^* 

The effective dielectric function ecff is obtained by com- 
parison of the reflectance coefficient of a surface of the 
composite materials under study to Fresnel's equations 
for equivalent homogeneous media.— We have found that 
a single value of ees can reproduce the angular depen- 
dence of the reflectance coefficients for all light polariza- 
tions with excellent numerical precision as long as the 
lattice period is small compared to the wavelength in the 
host material. 
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III. RESULTS AND DISCUSSION 

We have applied the above formahsm to study the 
following systems: periodic arrays of aluminum spheres 
placed in vacuum in simple cubic and fee configurations 
(Fig. [3]) , and fee arrays of silver spheres surrounded by 
either vacuum or silicon (Fig. [?]). The dielectric con- 
stants of the materials under consideration have been 
taken from optical data4^ In all cases, the distance be- 
tween nearest neighbor sites has been taken as 6 nm, 
although our results are quite insensitive to the choice 
of this parameter, as long as it is much smaller than the 
wavelength. 

Figures [3] and S] show contour plots of the so-called 
loss function, Im{— 1/eoff}, illustrating its dependence 
on metal filling fraction / and light frequency. In the 
absence of absorption, the loss function diverges at the 
plasma frequencies (ecff — 0), but in the actual, lossy 
systems under discussion this divergence is turned into 
finite peaks. Therefore, the regions of large values of the 
loss function (bright areas) correspond to plasmon exci- 
tations of the composites. Let us mention that the loss 
function inherits its name from the relevant role that is 
plays in electron energy loss spectroscopy (EELS) , where 
it provides the loss spectrum for electrons traversing the 
bulk of a material 4^ 

Aluminum is a good example of a material where the 
Drude formula ([T]) describes extremely well its bulk di- 
electric function with the parameters given in the cap- 
tion of Fig. [51 Therefore, we consider first a system of Al 
spheres in vacuum as a prototype system to study plas- 
mon chemistry in metallic particle arrays. We observe 
that there are more than one plasmon for each value of 
the filling fraction, unlike the single plasmon prediction 
of Eq. ^ based upon MGT. The plasmon frequency de- 
rived from that formula (solid curves in Fig. [3]) follows 
relatively well the highest-frequency plasmon of the rigor- 
ous calculation. This agreement is better in the fee lattice 
[Fig. [3Jb)]. However, this structure exhibit other plas- 
mons that originate in the splitting and mixing of indi- 
vidual particle modes. These plasmon frequencies merge 
into the isolated particle limit uip/ \/3 as / — > 0. For small 
values of /, there are only two plasmon branches that 
contribute significantly, arising from dipole-dipole inter- 
action among the spheres, although the low-frequency 
branch is less clear in the simple cubic lattice. For larger 
/, new plasmon modes show up in the loss function as 
a result of interaction between multiples of higher order, 
in agreement with Ref. This trend becomes dramatic 
near the percolation limit, where the material is expected 
to acquire good metallic behavior 

For an fee array of Ag nanoparticles (Fig. 2]), one can 
draw similar conclusions. When the host material is Si 
[Fig.Ufb)] the plasmon peaks are red shifted as compared 
to host vacuum [Fig. HJ^a)] and the plasmon peaks are 
relatively more spaced in frequency for a given value of /. 
Metallic nanoparticles embedded in Si constitute systems 
that have become recently realizable in practice,'^* so that 
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FIG. 5: Real and imaginary parts of the effective dielectric 
function tcff for an fee arrangement of Ag spheres surrounded 
by Si (solid and broken curves, respectively). Two different 
filling fractions of the metal have been considered: (a) / = 0.2 
and (b) / = 0.5. The dielectric function has been divided by 
50 in the low-u region to improve readability. 



the results presented here bear predictions that can be 
addressed experimentally. 

Both real and imaginary parts of the dielectric func- 
tion have been represented in Fig. [5] for fee lattices of Ag 
spheres in Si and with different values of the metal filling 
fraction (/ = 0.2 and / = 0.5). This figure indicates that 
the highest-energy plasmon branch observed in Fig. Iljb) 
arises from to the vanishing of the real part of eefi, with 
no further features showing up in the imaginary part, 
which is consistent with the fact that this plasmon is rel- 
atively close to the single plasmon predicted by MGT. 
This is quite different from the lowest-energy plasmon, 
where Ceff itself exhibits a resonant behavior translated 
into a Lorentzian-type shape, including a peak in the 
imaginary part. Other resonances of multipolar origin at 
intermediate energies [see Fig. \Mjo)] are mainly associ- 
ated to peaks in the imaginary part, which are particu- 
larly clear in Fig. Eljb) for / — 0.5. Our calculations for 
the rest of the combinations of lattice, metal, and embed- 
ding material considered in this work allow us to draw a 
similar picture in those cases. 

In conclusion, we have shown that the plasmons ex- 
hibited by metallic nanoparticles do hybridize when the 
particles are close to each other, and this produces new 
plasmon modes whose frequency can be tuned by chang- 
ing the distance between the particles. This suggests the 
possibility of designing temperature and pressure sensors 
based upon changes in the plasmon frequency with small 
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variations of the filling fraction, particularly near the per- 
colation limit (e.g., using short thiols to coat the parti- 
cles in order to have small separations between them). 
Our study demonstrates the failure of Maxwell-Garnett 
theory,^' that predicts only a single plasmon in these 
structures, rather than the rich structure of plasmons 
revealed by our rigorous solution of the full electromag- 
netic problem, not to mention Bruggeman's theoryf^iii^ 
that does not predict a plasmon at all.'*- 



sions. This work was supported by the University of the 
Basque Country UPV/EHU (contract No. 00206.215- 
13639/2001), by the late Spanish Ministerio de Cien- 
cia y Tecnologfa and the new Ministerio de Educacion 
y Ciencia (contract Nos. MAT2001-0946 and FIS2004- 
06490-C03-02, respectively), and by the European Union 
{Metamorphose Network of Excellence, NMP3-CT-2004- 
500252). SR was partially supported by the Emil Aalto- 
nen Foundation. 



Acknowledgments 

FJGA would like to thank F. Meseguer, M. KaU, 
J. M. Pitarke, and J. Fernandez for helpful discus- 



Corresponding author. 



^ R. A. Shelby, D. R. Smith, and S. Schultz, Science 292, 77 
(2001). 

2 J. B. Pendry, Phys. Rev. Lett. 85, 3966 (2000). 
^ T. Ung, L. M. Liz-Marzan, and P. Mulvaney, J. Phys. 
Chem. 105, 3441 (2001). 

* F. Caruso, M. Spasova, V. Salgueirio-Maceira, and L. M. 
Liz-Marzan, Adv. Mater. 13, 1090 (2001); L. M. Liz- 
Marzan and P. Mulvaney, J. Phys. Chem. B 107, 7312 

(2003) . 

^ V. F. Puntes, K. M. Krishnan, and A. P. Alivisatos, Science 

291, 2115 (2001). 
^ Y. Ding, Y. Kim, and J. Erlebacher, Adv. Mater. 16, 1897 

(2004) . 

J. C. Maxwell, Electricity and Magnetism (Oxford, Claren- 
don, 1873). 

* J. C. Maxwell-Garnett, Philos. Trans. R. Soc. London A 
203, 385 (1904); 205, 237 (1906). 

^ D. A. G. Bruggeman, Ann. Phys. (Leipzig) 24, 636 (1935). 
^° D. J. Bergman, Phys. Rev. B 19, 2359 (1979). 
" P. Sheng, Phys. Rev. Lett. 45, 60 (1980). 

W. Lamb, D. M. Wood, and N. W. Ashcroft, Phys. Rev. 

B 21, 2248 (1980). 

G. W. Milton, Appl. Phys. A 26, 1207 (1981). 

R. Tao, Z. Chen, and P. Sheng, Phys. Rev. B 41, 2417 

(1990). 

" D. J. Bergman and K.-J. Dunn, Phys. Rev. B 45, 13262 
(1992). 

" S. Datta, C. T. Chan, K. M. Ho, and C. M. Soukoulis, 
Phys. Rev. B 48, 14936 (1993). 

L. Martm-Moreno and J. B. Pendry, Nucl. lustrum. Meth- 
ods Phys. Res. B 96, 565 (1995). 

W. Lii, J. Dong, and Z.-Y. Li, Phys. Rev. B 63, 33401 
(2000). 

Z. Liu, C. T. Chan, P. Sheng, A. L. Goertzen, and J. H. 
Page, Phys. Rev. B 62, 2446 (2000). 
^° H. Ma, B. Zhang, W. Y. Tam, and P. Sheng, Phys. Rev. 
B 61, 962 (2000). 



J. E. Inglesfield, J. M. Pitarke, and R. Kemp, Phys. Rev. 
B 69, 233103 (2004). 

F. J. Garcia de Abajo and L. A. Blanco, Phys. Rev. B 67, 
125108 (2003). 

R. L Cukier, J. Karkheck, S. Kumar, and S. Y. Sheu, Phys. 
Rev. B 41, 1630 (1990). 

S.-Y. Sheu, S. Kumar, and R. L Cukier, Phys. Rev. B 42, 
1431 (1990). 

25 K. Hinsen and B. U. Felderhof, Phys. Rev. B 46, 12955 
(1992). 

R. G. Barrera and R. Fuchs, Phys. Rev. B 52, 3256 (1995). 
2^ C. Noguez and R. G. Barrera, Phys. Rev. B 57, 302 (1998). 
2® G. Mie, Ann. Phys. 25, 377 (1908). 

2^ J. R. Krenn, G. Schider, W. Rechberger, B. Lamprecht, 
A. Leitner, and F. R. Ausseneg, Appl. Phys. Lett. 77, 3379 
(2000). 

^" J. Aizpurua, P. Hanarp, D. S. Sutherland, M. KaU, G. W. 
Bryant, and F. J. Garcia de Abajo, Phys. Rev. Lett. 90, 
057401 (2003). 

E. Prodan, C. Radloff, N. J. Halas, and P. Nordlander, 
Science 302, 419 (2003). 
^2 P. Nordlander, C. Oubre, E. Prodan, K. Li, and M. L 
Stockman, Nano Lett. 4, 899 (2004). 

A. Taleb, V. Pussier, A. Courty, and M. P. Pileni, Phys. 

Rev. B 59, 13350 (1999). 
^'^ J. B. Pendry, A. J. Holden, W. J. Stewart, and L Youngs, 

Phys. Rev. Lett. 76, 4773 (1996). 
^5 V. C. Sundar, H. J. Eisler, and M. G. Bawendi, Adv. 

Mater. 14, 739 (2002). 

M. Gadenne, V. Podolskiy, P. Gadenne, P. Sheng, and 
V. M. Shalaev, Europhys. Lett. 53, 364 (2001). 
D. Zanchet, M. S. Moreno, and D. Ugarte, Phys. Rev. Lett. 
82, 5277 (1999). 

H. Fan, K. Yang, D. M. Boye, T. Sigmon, K. J. Malloy, 
H. Xu, G. P. Lopez, and C. J. Brinker, Science 304, 567 
(2004). 

^® Here, we work in the electrostatic limit, but the plasma fre- 



6 



quencies of spheres comparable in size to the wavelength 
are gcneraly red-shifted, as demonstrated via electron- 
induced light emission experiments using an electron mi- 
croscope in N. Yamamoto, K. Araya, and F. J. Garcia de 
Abajo, Phys. Rev. B 64, 205419 (2001). 
F. J. Garcia de Abajo, Phys. Rev. Lett. 82, 2776 (1999); 
Phys. Rev. B 60, 6086 (1999). 

F. E. Low, Classical Field Theory: Electromagnetism and 
Gravitation (John Wiley & Sons, New York, 1997). 
J. D. Jackson, Classical Electrodynamics (Wiley, New 
York, 1975). 

N. Stefanou, V. Yannopapas, and A. Modinos, Comput. 

Phys. Commun. 113, 49 (1998); 132, 189 (2000). 

A. Pack, M. Hietschold, and R. Wannemacher, Opt. 



Comm. 194, 277 (2001). 

E. D. Palik, Handbook of Optical Constants of Solids (Aca- 
demic Press, New York, 1985). 

D. W. McComb and A. Howie, Ultramicroscopy. 34, 84 
(1990); Nucl. Instrum. Methods Phys. Res. B 96, 569 
(1995). 

C. Pecharroman and J. S. Moya, Adv. Mater. 12, 294 

(2000) . 

C. Pecharroman, F. Esteban-Betegon, J. F. Bartolome, 
S. Lopez-Esteban, and J. S. Moya, Adv. Mater. 13, 1541 

(2001) . 

Bruggeman's theory leads to an inequality when one im- 
poses eefF = 0. 



